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We use Density Functional Theory to study interacting spinless electrons on a one-dimensional 
. . . quantum ring in the density range where the system undergoes Wigner crystallization. The Wigner 

C^^ ' transition leads to a drastic "collective" electron localization due to the Wigner crystal pinning, 

^3 \ provided a weak impurity potential is applied. To reveal this localization we examine a persistent 

^^ ■ current in a ring penetrated by a magnetic flux. Using the DFT-OEP method we calculated the 

Cn ' current as a function of the interaction parameter rs . We find that in the limit of vanishing impurity 

^ ' potential the persistent current stays constant up to a critical value of rg ~ 2.05 but shows a drastic 

Q , exponential decay for larger rs which reflects a formation of a pinned Wigner crystal. Above rg the 

^^ ' amplitude of the electron density oscillations exactly follows the (rs — rg) ' ^ behaviour, confirming 

a second-order phase transition as expected in the mean-field-type OEP approximation. 
0\ 

PACS numbers: 73.21.-b, 73.23.Ra 



CD 



cd 



O 

o 



> 

en 



1^ 

O 



X 



I. INTRODUCTION 

In the last years, the fabrication of quasi-onc- 
dimensional quantum rings became possibleii^. In such 
systems only few transverse states are occupied and by 
increasing the curvature of the confining potential the 
system can be made effectively one-dimensional. The 
number of electrons on the ring can be controlled by the 
gate electrode. The experimental studies of the rings 
with only one or two electrons were reported by Lorke et 
al.^ The possibility to vary the number of particles from 
very few to several hundreds enables experimentalists to 
tune the electron-electron interaction in a wide range. 
One of the most striking consequences of the interaction 
is the formation of a Wigner crystal^, a many-body state 
with electrons localized at discrete lattice sites. Yet it 
is well known that in an infinite one-dimensional sys- 
tem the fluctuations destroy the long-range order—. This 
raised doubts about the existence of a one-dimensional 
Wigner crystal which has become a long-debated sub- 
ject. Only in the nineties Glazman et al^ have shown 
that the arbitrarily weak pinning potential stabilizes the 
one-dimensional Wigner crystal. It was proven that the 
pinning potential suppresses the long-wavelength fluctua- 
tion modes which are responsible for destroying the long- 
range order. Due to the pinning potential the Wigner 
state is always localized in contrast to the electron liq- 
uid state^. Thus in the presence of a weak impurity the 
Wigner transition should manifest itself as electron local- 
ization. 

The critical r^ for a ID system estimated in the work 
of Glazman et al. was of the order of unity. For 
the two-dimensional electron gas Tanatar and Ceper- 
ley found a critical value rg = 37 ± 5, using a Monte- 
Carlo technique^. The reason for this large value is a 
very small shear modulus of the two-dimensional Wigner 
crystal^. In three dimensions a Wigner crystal is ex- 
pected at rs > 65 ± 10 (Ref. i). 

Electron localization seems to be a convenient signa- 



ture to observe a formation of the pinned one-dimensional 
Wigner crystal. However, in numerical simulations it is 
not quite evident how to quantify the localization of a 
correlated many body state. Several indirect criteria such 
as the inverse participation number— or the curvature of 
the ground state energjsi^ have been suggested to distin- 
guish between a localized and a delocalized state. How- 
ever, to the best of our knowledge, the electrons' ability 
to carry electric current - which is the most direct in- 
dication of the delocalized vs. localized behaviour - has 
not yet been explored. In this work we calculate the per- 
sistent current in a one-dimensional quantum ring pen- 
etrated by a magnetic flux. We apply a weak impurity 
potential which pins a Wigner state but practically does 
not influence the electron liquid state. 

In the density range where a Wigner crystal already 
exists as a ground state the persistent current has been 
studied analytically by Krive et ali^ for smooth po- 
tentials allowing semiclassical treatment. In a perfect 
ring the Wigner crystal rotates as a whole producing ex- 
actly the same current as non-interacting electrons. In 
the presence of a weak impurity potential the persistent 
current was found to be suppressed exponentially with 
the increasing impurity strength or the Wigner crystal 
stiffness ^-'^. 

We use Density Functional Theory (DFT) to calcu- 
late self-consistently the persistent current in a one- 
dimensional system with ten electrons. In the limit of 
vanishing (on a scale of the inter-electron Coulomb re- 
pulsion) repulsive potential we find that the current is 
independent of rg for rg < 2.05. At larger rg the persis- 
tent current decreases exponentially with increasing rg, 
indicating a localization of the electrons. At the tran- 
sition point the system undergoes a second-order phase 
transition which can be seen by considering the ampli- 
tude of the density oscillations S as an order parameter. 
We find that in a crystalline phase S exactly follows a 

1/2 

square root behaviour 6 ^ (rg — ?"§) . The stronger 
pinning potentials smear the phase transition such that 



no distinct transition point can be observed. 

The article is organized as follows. In section II we 
introduce the model of a one-dimensional quantum ring 
with the Gaussian impurity potential. We briefly dis- 
cuss the OEP approximationi^ii^ which is used for the 
exchange potential. We also introduce an Electron Lo- 
calization Function^* which is helpful for a real-space vi- 
sualization of the electron localization. In section III 
we describe the computational method for solving the 
self-consistent Kohn-Sham equations. In section IV we 
present our results for the persistent current as a function 
of rs . We consider impurity potentials of various ampli- 
tude and width and show how these parameters influence 
the current. The conclusions are given in section V. 



II. THE MODEL 

We study a system of A^ = 10 interacting spinless 
electrons in a one dimensional ring of circumference 
L = 2'kR. The ring geometry is accounted for via pe- 
riodic boundary conditions and x = LpR denotes the co- 
ordinate along the ring. A persistent current is induced 
by a vector potential A = {Ar,A^p) with a tangential 
component 



common choice is the "exact exchange" functional 
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that provides a magnetic flux $ through the ring. The 
vector potential is chosen such that the electrons move 
in a field- free space. 

Additionally, we introduce a repulsive Gaussian poten- 
tial centered at Xq 



Mmp(a;) = Vbexp 
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which should pin the Wigner crystal phase. 

We calculate the ground state current density for a 
given value of the magnetic flux and for a given strength 
and width of the impurity potential using Density Func- 
tional Theory. The self consistent Kohn-Shanti^ equa- 
tions for this system are given by 
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where index i labels the Kohn-Sham orbitals ipi and 
the eigenvalues e^. The electron-electron interaction is 
described by an effective one-particle scalar potential 
Mnt = Vh + V'oEP- Here, Vh is the Hartree potential 
and Vq^p is the exchange contribution. The latter is cal- 
culated in the KLI versioni^iii of the OEP methodi^iH. 
The central assumption of the OEP method is that 
the exchange-correlation energy functional can be writ- 
ten explicitely in terms of the Kohn-Sham orbitals. A 
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which has the form of the Fock energy but the wave- 
functions ^pi are the Kohn-Sham orbitals rather than the 
Hartree-Fock orbitals. Minimization of the full energy 
functional with respect to the density leads to an inte- 
gral equation for the exchange-correlation potential. In 
this work we use the exact-exchange functional and ap- 
ply the KLI approximation which allows to transform the 
OEP integral equation into a considerably simpler alge- 
braic equation. Still, it retains important features of the 
exact xc potential such as the derivative discontinuities 
and correct asymptotic behaviour—. 

Since DFT in the Kohn-Sham formulation is essentially 
a mean-field theory, fluctuations are not accounted for in 
our calculations. It is well known that fluctuations are 
particularly important in one dimension'^ . But since even 
an infinite one-dimensional Wigner crystal is stabilized 
by an arbitrarily weak pinning potential^ we expect that 
the fluctuations are effectively suppressed not only due 
to the pinning potential, but also due to the finite size of 
the ring. 

Whether the ground state of a many electron system is 
an electron gas-like one or a Wigner crystal state depends 
on the ratio of the kinetic energy and the Coulomb en- 
ergy. In one dimension this ratio is simply proportional 
to the electron density n, whereas in two and three di- 
mensions it is proportional to \pn and n '^ , respectively. 
Hence for high densities the kinetic energy dominates and 
the ground state is electron gas-like whereas for low den- 
sities the Coulomb repulsion favours the crystalline state. 

Experimentally it is most straightforward to vary the 
electron density to switch between weakly and strongly 
interacting regimes. Yet the variation of the electron 
number should alter the persistent current even in a non- 
interacting system which conceals the interaction effects. 
As we use a "persistent current criterion" to identify the 
Wigner transition we prefer to exclude the aforemen- 
tioned trivial single particle contribution and to retain 
only the influence of many-body effects. It can be done 
using an alternative (though somewhat artificial) way 
of controlling the ratio of kinetic and Coulomb energy. 
Namely, let us consider the effective electron mass m* as 
a free parameter. In one dimension, the energy ratio rg 
is proportional to m*: 
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where ob and m^ are the Bohr radius and the "true" ef- 
fective electron mass in the host material. The persistent 



current density 
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should be calculated with the fixed "true" effective elec- 
tron mass TOq. Here, $o = - is the flux quantum and 
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is the density. 

Also, the ratio of the kinetic energy to the impurity 
potential must be kept constant when changing rs via 
changing m* . Otherwise the current density of a system 
of non- interacting electrons would depend on r^- The 
impurity potential strength Vq must be renormalized as 



Vo ^ Vo* = Vo- 



(8) 



The potential renormalization ([5]) guarantees that the 
artificial variation of the electron mass results in a de- 
pendence of the persistent current on rg solely due to 
the electron-electron interaction. 

Equation ^ expresses the current density via the 
Kohn-Sham orbitals within the framework of the ordi- 
nary density-based DFT. It is not, however, strictly jus- 
tified since the common DFT Kohn-Sham equations by 
construction yield the exact ground state density but not 
the current density. Strictly speaking, one has to employ 
the current density functional theory^^ (CDFT) which 
expresses the ground-state energy functional as a func- 
tional of the density and the paramagnetic current den- 
sity. The Kohn-Sham orbitals in CDFT thus give the ex- 
act current density of the interacting system. However, 
the CDFT corrections are, as a matter of fact, usually 
very small. For example in a recent paper^*^ is shown 
that the orbital magnetic moments in magnetic (Fe, Co 
and Ni) and non-magnetic (Si and Ge) solids calculated 
with CDFT only slightly differ from those calculated with 
DFT. Hence we expect that Eq. (jG)) evaluates the current 
reasonably well, keeping in mind, that for our purposes 
not the current value itself, but its critical rs-dependence 
close to the Wigner transition is of interest. 

In addition to the current density we also use the Elec- 
tron Localization Function (ELF) ^'^ to visualize the elec- 
trons' localization. The idea behind the definition of 
the ELF is that the more localized electron produces a 
stronger repulsion of the other like-spin electrons due to 
the Pauli exclusion principle. According to this picture 
the ELF measures the probability to find a second elec- 
tron (with the parallel spin) anywhere close to a refer- 
ence electron. It is defined such that its value of one half 
means a homogeneous electron-gas like state whereas a 



value of one refers to a perfectly localized electron at this 
point in space. 

In its original dcfinitionii the ELF was formulated for 
the real wavefunctions only. Recently it was generalized 
to the time-dependent case^^ where complex wavefunc- 
tions have to be employed. This form of the ELF is also 
suitable for the current-carrying static system we con- 
sider. It is given by 
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In this expression t{x) — ^ '^^ \dx'fi{x)\'^ is the kinetic 
energy density of the Kohn-Sham system and t'^°™(x) = 
g ^, n'^(x) is the respective quantity in a one-dimensional 
homogeneous electron gas with density n{x). 



III. COMPUTATIONAL METHOD 

For numerical solution of the Kohn-Sham equations 
^ we use a real space method. We expand the wave 
functions ipi{x) using a spline basis^^ 
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with the complex coefficients a^ and the real basis func- 
tions 
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(12) 
The spline nodes are Xy and h is the distance between 
the two adjacent nodes. The basis functions (|12p are not 
orthogonal which means that the overlap matrix 



S^,.v = / dxbf^{x)by{x) 



(13) 



is not diagonal. With this representation of the wave 
functions, the Schrodinger equation reads 
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with the Hamiltonian matrix 

Hf,,,= [ dxb^ix)Hb,{x). (15) 



At the first step this generahzed eigenvalue equation is 
transformed into a standard eigenvalue equation. We use 
a Cholesky decomposition^^ of the overlap matrix 
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into a lower triangular matrix L and its transpose and 
write the eigenvalue equation as 
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The matrix L ^H{l^\ is diagonalized using the 

zheev-routine form the LAPACK library^"* and the re- 
sulting eigenvector L^ai is transformed back to obtain 
the eigenvector Si of the original generalized eigenvalue 
problem. 

The starting point for the iterative self-consistent pro- 
cedure is a system of non-interacting particles i.e. a sys- 
tem with Vh = V^ = 0. The resulting non- interacting 
eigenfunctions are then used to construct the first ap- 
proximation for the Hartee- and the exchange potential. 
In the subsequent iterations the Hartree- and the ex- 
change potential are calculated from the eigenfunctions 
of the previous step^. As a measure of the convergence 
we consider the maximum difference between two Kohn- 
Sham eigenvalues in the n-th and (n — l)-th iteration 
step: 
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We found that this difference has to be extremely small 
compared to the Kohn-Sham eigenvalues themselves 
which are of the order of several tens of meV, namely 
A « 10~"'^°meV. The reason for this very small num- 
ber are low energy excitations which correspond to a 
charge displacement over a large distance in the system. 
If the chosen A is too large, one encounters a density 
range where the system seems to be in a delocalized state 
whereas in fact it becomes localized after the solution is 
converged. Generally, a very high computational accu- 
racy is required to distinguish correctly between a local- 
ized and a delocalized state of the system. 



IV. RESULTS 

In this section we present the results of our calculations 
of the persistent current in the one-dimensional quantum 
ring. For the effective electron mass and the dielectric 
constant we have chosen the GaAs values jtiq — 0.0665toc 
and e = 12.5. The value of the magnetic field flux was 
chosen as $ = 0.3$o. In fact, the particular magnitude of 
the flux does not matter provided the current distinctly 
exceeds numerical inaccuracy. 

For the Wigner crystal pinning we apply a narrow im- 
purity potential of a width a — 0.025L much smaller 
than the average distance between electrons jr = O.IL. 
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FIG. 1: (Color online) The persistent current as a function 
of rs for a Gaussian impurity potential with a half maximum 
width of 2.5% of the ring circumference. The current is nor- 
malized to its value jo in a non-interacting system and poten- 
tial strength Vo = 10~^meV. The long-dashed line -'/jq — 1 
corresponds to the interaction-free system. 
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FIG. 2: (Color online) Logarithmic version of the plot in Fig. 
[T] The exponential dependence of the persistent current on 
rs is clearly seen. 



The persistent current is calculated as a function of rs, 
the latter being altered by varying m*, according to Eqs. 
(O, ([8]). The current is normalized to its value jo for 
non-interacting electrons in the presence of an impurity 
potential with unrenormalized strength Vq = lO^^meV. 
The results for various impurity potential strengths are 
shown in Fig. fTl The dashed line -^ — 1 reflects the 

Jo ^ 

current independence of rg for noninteracting electrons. 

As seen in Fig. [T]for the smallest Vb — lO^'^ meV, one 
can clearly distinguish two different regions of rg . Below 
the critical value of rg « 2.05, the persistent current is 
independent of rg. Its magnitude is the same as in the 
non-interacting system which means that the interacting 
system is electron gas-like. In contrast, for rg > rg, the 
persistent current drops exponentially with increasing rg 
which is seen explicitely from the linear dependence of 
log-* /jo on rg shown in Fig. [21 This signifies the forma- 
tion of the Wigner crystal pinned by an extremely small 
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FIG. 3: (Color online) Electron Localization Function in the 
presence of a weak (Vb = O.OOlmeV) potential. Shown is 
the ELF for different values of rs. Solid line: rs — 0.1, 
dashed line: rs — 2.06, dotted line: rs — 2.5, dash-dotted 
line: rs = 5.0. An ELF value of one corresponds to perfect 
localization whereas an ELF value of one half means homo- 
geneous electron gas-like delocalization. 



impurity potentiaL Hence the value r^ = 2.05 can be 
interpreted as a critical rs of the Wigner transition. 

This interpretation is supported by the ELF plot in 
Fig. [3] For rs < 2.05 we find an ELF value of one 
half, corresponding to completely delocalized electrons. 
This changes drastically when rs exceeds r^. With in- 
creasing rs the electrons tend to localize at discrete lat- 
tice sites. At rs « 5 they arrange in an "almost classi- 
cal" one-dimensional lattice. The complete localization 
is achieved within a rather narrow interval of rs as ex- 
emplified in Fig. [3] by the ELF graphs for rs = 2.06 
and rs = 2.5. This refiects the exponential decay of the 
persistent current as shown in Fig. [T] 

We believe, that within our numerical accuracy the 
solid curve in Fig. [1] corresponds to the case of the 
"vanishing" external potential. Such a potential does 
not disturb the Wigner transition, but provides the pin- 
ning. The particular potential strength and width should 
be then unimportant. We tested this calculating the 
current density for several values of the width of the 
pinning potential (all with Vq = 10~'^meV) and found 
that the persistent current follows exactly the same rs- 
dependence. However, the convergence is getting much 
harder for wider potentials since the "smoother" poten- 
tials are less effective in pinning the Wigner crystal. For 
Vq values below lO^'^mcV the convergence could not be 
reached. Yet using a semiclassical approach— it can be 
shown analytically that the current value jo of a non- 



interacting system is indeed recovered for T^n 



0. 



The critical rg — 2.05 we obtained in this work is 
of the same order as the values for rg found in a pre- 
vious work2^ for a different model using the ground 
state energy curvature^ as a localization criterion. In 
the presence of a disorder potential with an amplitude 
Ay = 0.02meV a Wigner transition has been observed 
in the range 2.08 < rg < 5.04 depending on the model 



FIG. 4: (Color online) Electron Localization Function in the 
presence of an intermediate (Vb = 5.0meV) narrow poten- 
tial. Shown is the ELF for different values of rs. Solid line: 
rs = 0.1, dashed line: rs = 1.5, dotted line: rs — 2.0, dash- 
dotted line: rs — 2.5, long-dashed line: rs = 5.0. For inter- 
mediate values of rs the electrons next to the impurity are 
more localized. This localization increases gradually with in- 
creasing rs. For large values of rs the ELF is the same as 
found in the case of the weak potential (see Fig [3]). 



for the electron-electron interaction. 

The other three curves in Fig. [T]show the current of the 
interacting system for Vb = l.OmeV, Vb — 5.0meV and 
Vb = lO.OmeV. Although at Vb = l.OmeV there is still 
the range of rs where j = jo, the sharp kink at rs = r^ 
vanishes. The transition smoothing is more pronounced 
for Vb — 5.0meV and Vb = lO.OmeV where no region of 
rs where the current is independent of rs is seen. 

It should be emphasized that the dependence of the 
normalized current on rs is solely due to the electron- 
electron interaction. The smooth decrease of the cur- 
rent with increasingly strong Coulomb interaction ob- 
served for stronger impurity potentials (Vq = 5.0meV 
and Vb = lO.OmeV) reflects a gradual localization of the 
many-body state instead of a distinct phase transition. 
This behaviour parallels the absence of a sharp phase 
transition in an external potential field that lowers the 
symmetry of the high-symmetry phased. 

An estimate of the Coulomb energy of two electrons at 
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shows that it is indeed of the order of the pinning poten- 
tial which smoothes out the phase transition and induces 
a gradual localization. For Vb > ImeV and at interme- 
diate values of rs it can be seen directly from the ELF 
plots (Fig. 3]) that the localization is more pronounced 
next to the pinning potential. This indicates that a grad- 
ual localization seen in Fig. [l]is driven by the interplay 
between the long-range Coulomb repulsion and the in- 
teraction with the short-range impurity potential, both 
being of the same order. 

The sharp transition we found for a "vanishing" im- 




FIG. 5: (Color online) Amplitude of the density oscillations 
as a function of rs for a weak impurity potential (Vb = 
O.OOlmeV). The solid black curve shows the calculated data, 

the dashed red curve is a square root (rs — r§) " behaviour. 



purity potential (solid line in Fig. [T]) is a second or- 
der phase transition from an electron liquid state to the 
Wigner crystal state. This can be verified by plotting the 

rs-dependence of the order parameter S which shows a 

1 / 
behaviour S ^ (^s — ^s) ^^ '^S > '"§' i-^- ^^ ^^^ ^O'^" 
symmetry phased. Indeed, taking the amplitude of the 
density oscillations as the order parameter S, we obtain 
an exact square root dependence, as shown in Fig. [S] 
The second-order type of the transition we observe in 
our calculations is quite natural for the mean-field-type 
DFT-OEP approach. 

From the exponential dependence of the current on rg 
(Fig. [1]) we can deduce the relation between the persis- 
tent current density and the order parameter 



j{^) = Joexp(-a(5^) 
where the numerical factor a = 0.033i^. 



(20) 



CONCLUSIONS 



In this article we investigated numerically the influence 
of the electron-electron interaction on the ground state 
of a one-dimensional electron gas confined in a ring ge- 
ometry. To break the rotational invariance of the ring 
we introduce a weak "impurity" potential. This poten- 
tial does not affect the delocalized electron liquid phase, 
but provides a pinning of the crystalline Wigner phase. 
We employ a persistent current in the ring as a measure 
of the Wigner crystal pinning. For a sufficiently weak 
impurity potential we found that for rs < r^ the cur- 
rent density of the interacting system is exactly the same 
as the current density of a non-interacting electron gas. 
For rs > r^ the current of the interacting system de- 
cays exponentially with increasing rg while the current 
of a non- interacting system remains constant. This be- 
haviour clearly shows the formation of the Wigner crys- 
tal in a one-dimensional system. This interpretation is 
confirmed by the ELF plots which reveal the delocalized 
electron distribution below the critical rg and a local- 
ized one above Tg. At rg = rg the system undergoes 
a second-order phase transition from an electron liquid 
to a Wigner crystal. This is evident from the square 
root dependence of the amplitude of the density oscil- 
lations (taken as the order parameter) on rg above the 
critical value. Experimentally, this transition should be 
observable as a sharp decrease of the ring's magnetiza- 
tion when the electron density is lowered. However, in a 
real experiment this transition will be superposed with 
the interaction-independent variation of the current den- 
sity due to the variation of the particle number. The 
critical value rg = 2.05 we find for the Wigner transition 
is consistent with the density rangc^ in which Glazman 
et al. expected the existence of a stable one-dimensional 
Wigner crystal. 
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